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ABSTRACT 


Impurities  such  as  H,  P,  S,  B,  etc.,  have  a  very  low  solubility  in  iron  and,  therefore, 
prefer  to  segregate  at  the  grain  boundaries  (GBs).  In  order  to  analyze  the  energetics 
of  impurities  on  the  iron  GB,  the  Linear  Muffin  Tin  Orbitals  (LMTO)  calculations  were 
performed  on  a  simple  8-atom  supercell  emulating  a  typical  (capped  trigonal  prism)  GB 
environment.  The  so-called  environment-sensitive  embedding  energies  (EE)  were  calcu¬ 
lated  for  H,  B,  C,  N,  O,  Al,  Si,  P,  and  S  as  a  function  of  the  electron  charge  density 
due  to  the  host  atoms  at  the  impurity  site.  It  was  shown  that,  at  the  electron  charge 
density  typical  of  a  GB,  B  and  C  have  the  lowest  energy  among  the  analyzed  impurities 
and,  thus,  would  compete  with  them  for  the  site  on  the  GB  tending  to  push  the  other 
impurities  off  the  GB.  The  above  energies  were  then  used  in  a  modified  Finnis-Sinclair 
embedded  atom  approach  for  calculating  the  equilibrium  interplanar  distances  in  the  vicin¬ 
ity  of  a  (111)  23  tilt  GB  plane,  both  for  the  clean  GB  and  that  with  an  impurity. 

These  distances  were  found  to  be  oscillating,  returning  to  the  equilibrium  spacing  be¬ 
tween  (111)  planes  in  bulk  BCC  iron  by  the  10th  to  12th  plane  off  the  GB  plane.  H, 

B,  C,  N,  and  O  actually  dampen  the  deformation  wave  (making  the  oscillation  ampli¬ 
tudes  less  than  in  the  clean  GB)  while  Al,  Si.  P,  and  S  result  in  an  increase  of  the 
oscillations.  The  effect  of  B,  C,  N,  and  O  may  be  interpreted  as  cohesion  enhance¬ 
ment;  this  conclusion  supports  our  earlier  first-principles  results  on  B  and  C. 
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INTRODUCTION 


The  reduced  cohesion  of  grain  boundaries  (GBs)  is  known  to  be  often  the  controlling 
factor  limiting  ductility  and,  actually,  the  performance  and  reliability  of  high-strength  metallic 
alloys.'  Intergranular  embrittlement  in  metals  is  usually  caused  by  impurities  segregating 
towards  the  GBs.^'^  Impurities  present  in  bulk  concentrations  of  10*^  to  10“^  atomic  percent 
can  result  in  a  dramatic  decrease  of  plasticity,  drastically  degrading  mechanical  properties  of 
metallic  alloys  and  thus  posing  significant  technological  problems.  This  detrimental  effect  of 
just  ppm  of  impurities  may  be  readily  understood:  a  simple  estimate  shows  that  a  ppm 
amount  of  impurity  is  sufficient  for  saturating  all  the  grain  boundaries  in  a  typical  grain-size 
polycrystal. 

Recent  progress  in  developing  efficient  methods  of  first-principles  calculations  and  compu¬ 
tational  algorithms  made  possible  systematic  studies  of  the  role  of  impurities  in  intergranular 
cohesion  on  the  electron-ion  level.  Calculations  on  both  cluster,^  two-dimensional,'’'’  and 
supercell^’'®  models  of  GBs  with  impurities  have  provided  an  in-depth  insight  into  mecha¬ 
nisms  of  GB  decohesion. 

Since  the  first-principles  electronic  calculations  on  low-symmetry  systems  (such  as  lattice 
defects  or  GBs)  are  still  extremely  complicated  and  expensive,  semi-empirical  methods  based 
upon  solid  first-principles  foundations  have  been  developed;  among  them,  the  most  popular  is 
the  Embedded  Atom  Method  (EAM).""*^  This  method  has  been  successfully  used  in  a  wide 
variety  of  calculations. 

The  purpose  of  this  report  is  to  elucidate  energetics  of  impurities  on  an  iron  GB  and 
analyze  the  effect  of  impurities  on  the  intergranular  cohesion  in  iron. 

Rather  than  doing  sophisticated  first-principles  calculations  on  multi-atom  models  of  a 
GB,  it  was  decided  to  calculate  the  quantity  which  may  be  called  environment-sensitive 
embedding  energy  (EE),  the  energy  of  an  impurity  atom  in  an  atomic  environment  typical 
for  a  GB  (capped  trigonal  prism).  Knowledge  of  these  energies  for  different  impurities 
enables  one  to  compare  the  relative  stability  of  a  particular  impurity  in  the  iron  GB.  The 
approach  taken  is  obviously  an  extension  of  the  so-called  Effective  Medium  (EM)  theory.'*'"^ 

Having  calculated  the  EEs  for  a  number  of  impurity  atoms,  one  can  use  this  information 
on  the  impurity  energetics  in  a  modified  EAM  approach  for  calculating  the  GB  relaxation. 
The  latter  calculations  enables  one  to  draw  important  conclusions  regarding  the  intergranular 
cohesion  in  iron  in  the  presence  of  a  definite  impurity  in  the  GB. 

RESULTS  AND  DISCUSSION 
Environment  Sensitive  Embedding  Energies 

The  basic  idea  of  the  EM  approach  was  to  replace  the  low  symmetry  system  consisting  of 
an  atom  plus  a  host  matrix  by  a  high  symmetry  effective  system  of  the  atom  and  the  homoge¬ 
neous  electron  gas  of  a  density  equal  to  that  seen  by  the  atom.  The  energy  of  interest  was 
called  the  embedding  energy  and  was  equal  to  the  energy  difference  between  the  atom  embed¬ 
ded  in  the  electron  gas  and,  separately,  the  isolated  atom  and  the  electron  gas. 
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The  EM  theory,  as  a  first  and  very  crude  approximation,  completely  neglected  any  cova¬ 
lent  effects,  though  one  could  expect  that  a  metalloid  impurity  in  a  transition  metal  would 
develop  strong  sp-d  hybridization  resulting  in  covalent  bonds.  The  introduction  of  covalent 
effects  via  perturbation  theory,**  resulted  in  significant  corrections  in  the  embedding  function 
for  hydrogen.  Further  attempts  to  improve  upon  the  EM  method  were  undertaken  in  recent 
years. 

Improvements  in  the  EM  theory  actually  made  calculations  more  sophisticated,  spoiling 
the  elegant  simplicity  of  the  original  method.  Rather  than  introducing  corrections  to  the  EM 
concept,  it  was  decided  to  perform  first-principles  calculations  on  a  simplified  model  of  a  GB 
environment,  varying  hydrostatically  all  the  characteristic  volumes,  thus  generating  a  series  of 
environment  sensitive  EEs  as  a  function  of  the  electron  charge  density  due  to  the  host  (iron) 
atom  at  the  impurity  site.  In  this  case  an  impurity  is  actually  embedded  into  a  crystal  lattice 
environment  rather  than  into  an  electron  jellium  as  in  the  original  EM  theory. 

The  model  chosen  for  the  GB  environment  is  an  8-atom  hexagonal  supercell  (Fe6X, 
where  X  is  an  impurity  atom).  The  supercell  is  shown  in  Figure  1  together  with  the  capped 
trigonal  prism  coordination  of  the  surrounding  iron  atoms. 


Figure  1.  The  FeeX  hexagonal  supercell  emulating  a  typical  trigonal 
prism  environment  of  iron  atoms  in  the  (1 1 1)  X3  GB:  a)  the  supercell; 
b)  the  trigonal  prism  coordination;  0  Fe.*  Impurity. 


A  trigonal  prism  GB  configuration  follows  from  the  theory  of  hard-sphere  packing.^’ 
Atomistic  relaxation  studies  have  also  shown^^  that  an  impurity  atom  such  as  P  or  B  is 
likely  to  occupy  an  interstitial  position  in  the  center  of  the  trigonal  prism  formed  by  iron 
atoms  in  the  GB  core  (even  if,  as  in  the  case  of  P  and  B,  the  impurity  forms  a  substitu¬ 
tional  solid  solution  with  the  host).  The  hexagonal  supcrcell  has  been  chosen  both  be¬ 
cause  of  its  relatively  high  symmetry  and  its  emulation  of  a  (111)  23  GB  environment.*^  ’** 

The  spin-polarized  scalar-relativistic  Linear  Muffin  Tin  Orbitals  (LMTO)"  calculations 
were  performed;  the  von  Barth-Hedin^"*  exchange-correlation  and  the  frozen  core 
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approximations  were  also  used.  First,  a  series  of  calculations  (for  six  different  volumes)  was 
performed  with  an  impurity  absent  from  the  supercell;  i.e.,  an  empty  sphere  of  the  same  ra¬ 
dius  as  that  of  the  radius  of  the  impurity’s  Wigner-Seitz  sphere  was  substituted  for  the  latter. 
Similar  calculations  were  then  performed  for  each  of  the  impurities;  H,  B,  C,  N,  O,  Al,  Si, 

P,  and  S.  The  EEs  were  defined  as  follows; 

EE  =  E(Fe6«)-E(Fe60)-E(»)  (1) 

where  E(Fe6«)  and  E(Fe60)  are,  respectively,  the  energies  of  the  supercell  with  and  without 
the  impurity  (O  stands  for  an  empty  sphere  substituted  for  the  impurity  atom),  and  E(«)  is 
the  energy  of  the  free  impurity  atom.  In  order  to  make  the  calculations  more  consistent,  the 
values  of  E(Fe6«)-E(Fe60)  extrapolated  to  the  zero  charge  density  (n  =  O),  were  chosen  to 
be  used  as  E(«)s  which  would  correspond  to  the  energies  of  impurities  in  the  GB  environ¬ 
ment  with  the  host  crystal  lattice  infinitely  expanded.  The  EE  energies  (see  Equation  1)  as  a 
function  of  n,  the  electron  charge  density  due  to  iron  atoms  at  the  impurity  site,  are  pre¬ 
sented  in  Figure  2. 


Figure  2.  The  environment  sensitive  EEs  versus  electron  charge  density  (in  atomic  units,  a.u.). 
a)  Hydrogen  and  impurities  of  the  second  period;  0)  Impurities  of  the  third  period. 


Plots  in  Figure  2  explain  an  experimentally  observed  phenomenon  known  as  site  competi¬ 
tion.  As  can  be  seen,  in  the  range  of  electron  charge  density  typical  of  a  GB  (0.015  a.u.  to 
0.025  a.u.),  B  and  C  have  the  lowest  energies  and,  thus,  would  push  the  other  impurities  off 
the  GB.  In  fact,  C  is  known  to  successfully  compete  with  P^^  and  and  N  -  with  P"^ 

and  B  should  even  be  more  efficient  in  pushing  the  harmful  impurities  off  the  GB. 

though  the  energy  difference  between  B  and  C  is  not  very  significant.  Since,  in  principle, 

EEs  are  to  be  sensitive  to  the  type  and  orientation  of  the  GB,  in  some  GBs  B  and  C  may 
compete  in  a  reverse  way.  Also,  it  can  be  seen  that  H  would  successfully  compete  with  all 
the  impurities  but  B  and  C;  N  being  a  marginal  case. 

The  plots  in  Figure  2  also  reveal  an  important  aspect  of  GB  impurity  behavior.  All  the 
plots  (except  for  H)  have  well  pronounced  minima.  The  positions  of  the  minima  correspond 
to  electron  density  at  the  impurity  site  due  to  the  surrounding  Fe  atoms  which  would  occur  if 
the  GB  were  allowed  to  relax  in  such  a  way  as  lo  minimize  the  impurity’s  energy.  The 
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minima  positions  systematically  shift  towards  lower  densities  with  the  impurity  loosing  its  com¬ 
petitive  power.  Lower  charge  density  means  a  more  loose  GB,  less  strong  and  more  prone  to 
decohesion.  The  minimization  of  the  total  GB  energy  (not  only  the  energy  of  the  impurity 
atom)  gives  the  characteristic  charge  densities  which  are  somewhat  higher  than  those  in  the 
minima.  However,  from  this  point  of  view,  only  S,  P,  Si,  and  AJ  are  the  obvious  candidates 
for  being  decohesive,  O  and  N  are  marginal,  while  B,  C,  and  H  may  be  called  cohesion 
enhancers  (the  behavior  of  H  will  be  discussed  in  more  detail  below).  Calculations  of  GB 
relaxation  shed  more  light  on  intergranular  cohesion/decohesion  due  to  impurities. 

Grain  Boundary  Relaxation 

First-principles  calculations  are  still  too  expensive  to  be  used  for  the  investigation  of 
GB  relaxation.  A  semi-empirical  EAM  is  obviously  more  appropriate.  In  order  to  find  both 
embedding  functions  and  pair  potentials,  the  EAM  uses  experimental  data  such  as  cohesive 
energies  and  elastic  moduli  for  the  system  of  interest.  This  approach  is  easily  applicable  to 
pure  metals  and  compounds  but  may  become,  in  fact,  inappropriate  if  the  effect  of  an  iso¬ 
lated  impurity  atom  or  an  impurity  atom  on  a  GB  is  to  be  studied.  For  example,  in  principle 
one  can,  in  a  metal-hydrogen  system,  use  the  embedding  function  and  the  pair  potentials 
found  from  experimental  information  on  that  metal’s  hydrides.  This  information,  however,  can 
be  misleading  since  the  effect  of  isolated  hydrogen  atoms  on  electronic  structure  and  cohesive 
properties  of  the  metal  may  be  completely  different  from  that  of  periodic  arrays  of  hydrogen 
atoms  typical  of  hydrides.  Krasko  and  Olson found  that  hydrogen  in  the  iron  GB  does  not 
contribute  its  electron  to  the  iron  d-valence  band  at  all,  contrary  to  a  general  belief  that  in 
transition  metals  the  hydrogen’s  electron  will  inevitably  go  to  a  d-band. 

In  order  to  resolve  this  difficulty  the  energy  contributions  due  to  impurity  atoms  in  the 
GB  were  calculated  by  using  the  EEs  discussed  above.  Since  the  EAM  functions  arc  also 
fundamentally  dependent  upon  the  electron  charge  density  at  an  atom  site,  the  EEs  may 
simply  be  added  to  the  EAM  energy  of  the  host  atoms: 

E  =  Sr  Eemb(n(R))  +  1/2Sr,r-  V(R,R')  EE(n(Rin,p))  (2) 

where  Eemb  (n)  and  V(R,R')  are  the  EAM  embedding  energy  and  the  pair  potential  as 
found  for  the  bulk  BCC  Fe  (the  Finnis-Sinclair  functions  and  parameters’^  were  used).  The 
third  term  is  the  energy  of  the  impurity  atom.  R  and  R'  are  the  positions  of  the  host 
atoms,  Rimp  is  that  of  the  impurity,  and  n(R)  and  n(Rjmp)  are  the  electron  charge  densities 
at  the  site  of  a  host  atom  and  the  impurity,  respectively.  The  electron  charge  density  at  a 
given  site  can  be  taken  to  be  a  superpusition  of  the  free  atom  charge  densities  or  found 
from  more  sophisticated  procedures. 

As  mentioned  above,  the  GB  environment  dealt  with  was  that  of  the  (111)  S3  tilt  GB. 

The  GB  structure  can  be  represented  as  a  succession  of  (111)  hexagonal  planes: 


....CBACBACBACBACBABCABCABCABCABC.... 


(the  GB  plane  is  marked  by  The  CBABC  atomic  structure  of  the  core  of  the  GB  (clean 
or  with  an  impurity)  was  just  emulated  by  the  supercell  and  shown  in  Figure  1.  In  order  to 
find  the  GB  structure  corresponding  to  a  minimum  of  energy,  (see  Equation  2)  the  inter- 
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An  interesting  feature  of  the  clean  GB  relaxation  is  that  the  distance  between  the  second 
and  third  planes  is  only  half  of  the  (111)  interplanar  distance  in  bulk  BCC  Fe  (0.407  A  versus 
0.825  A).  The  LMTO  calculations^^  show  that  it  is  only  the  interaction  between  Fe  atoms  in 
planes  3,  4,  5,  and  6  that  prevents  plane  3  from  collapsing  into  plane  2  (the  <y-phase  configura¬ 
tion).  The  site- projected  electronic  densities  of  states  of  the  Fe  atoms  in  planes  1  through  3  are 
very  similar  to  that  typical  of  an  <w-phase.*® 
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The  impurity  atoms,  H,  B,  C,  N,  and  O  result  in  some  damping  of  the  relaxation  deforma¬ 
tion  waves;  i.e.,  decreasing  the  oscillation  amplitudes.  This  damping  is  most  pronounced  for 
H,  B,  and  C  while  N  and  O  are  marginal.  The  amplitude  increases  with  Al,  Si,  P,  and  S. 
Damping  the  deformation  wave  may  be  interpreted  as  cohesion  enhancement  while  increasing 
the  deformation  wave  oscillations  may  be  thought  of  as  resulting  in  decohesion. 

C  and  N  are  known  to  be  GB  cohesion  enhancers  in  steels;  *  O,  Si,  P,  and  S,  being 
strong  embrittlers,  are  believed  to  weaken  the  GB  cohesion. The  results  agree  with  both 
this  experimental  information  and  the  earlier  first-principles  calculations^  on  B,  C,  P,  and  S. 

H  is  one  of  the  worst  embrittlers  of  steels.  At  the  same  time,  the  first-principles  calcula¬ 
tions*®  do  not  show  any  dramatic  GB  decohesion  due  to  H.  In  fact,  H  is  quite  neutral;  a 
weak  H-Fe  covalent  bond  does  not  significantly  disturb  the  GB  energetics,  and  the  H  electron 
does  not  participate  in  metallic  bonding.  Also,  the  calculations  of  both  the  EE  and  the  GB 
relaxation  do  not  show  any  sign  that  the  presence  of  H  in  the  GB  should  weaken  it  dramati¬ 
cally.  However,  the  EAM  calculation*^  has  shown  that  a  H  atom  on  a  Ni  GB  does  weaken 
the  metallic  bond  across  the  GB,  lowering  the  fracture  stress  by  some  15%. 

From  a  thermodynamic  point  of  view,  the  impurity’s  embrittling  potency  depends  upon  the 
difference  between  the  free  energies  of  the  impurity’s  segregation  on  the  initial  GB  and  on  the 
two  free  surfaces  emerging  upon  fracture.  The  higher  the  difference,  the  stronger  the  embrit¬ 
tling  potency  of  the  impurity.  As  a  less  rigorous  but  simpler  criterion,  in  Reference  30  the  sub- 
limination  energy  differences  between  the  host  and  impurity  were  calculated  in  an  ideal  solution 
model  for  over  60  elements.  The  effects  of  impurities  on  GB  stability  can  be  analyzed  by  simply 
comparing  the  GB  energy  differences  between  the  GB  with  impurities  and  the  clean  GB. 

The  corresponding  values  for  the  impurities  discussed  arc  plotted  in  Figure  4.  It  is  noted 
that  the  AE  values  are  negative  for  impurities  that  can  be  thought  of  as  cohesion  enhancers, 
and  positive  for  the  embrittlers.  O  and  Al,  the  marginal  impurities,  have  low  negative  values. 
The  AE  for  H  is  also  slightly  negative. 


Figure  4.  AE  =  EobM-EgbIcIh),  the  energy  difference  tjefween  the  GB  with 
impurity  X  and  clean  GB,  versus  the  Periodic  Chart  group  number. 


6 


CONCLUSIONS 


The  environment  sensitive  EEs  (the  energies  of  an  impurity  atom  in  an  environment  typi- 
:al  of  a  GB  in  Fe  as  a  function  of  electron  charge  density  at  the  impurity  site)  were  calcu- 
ated  for  a  number  of  impurities  often  present  in  Fe  and  Fe-base  alloys.  The  knowledge  of 
;hese  energies  enabled  us  to  explain  the  well  known  site  competition  effect  and  draw  import- 
mt  conclusions  about  the  relative  stability  of  GBs  with  different  impurities.  The  calculation 
Df  GB  relaxation  also  revealed  some  interesting  features  of  the  GB.  The  impurities  known  as 
:ohesion  enhancers  tend  to  suppress  the  deformation  wave  oscillations  while  the  embrittling 
impurities  increase  the  oscillations.  The  detrimental  effect  of  hydrogen  in  steels  cannot  be 
explained  from  the  knowledge  of  its  static  energetics.  In  order  to  explain  hydrogen’s  extreme 
jmbrittling  potency,  one  should  probably  investigate  the  dynamics  of  hydrogen  atoms  in  the 
GB  vicinity.  Since  the  EEs  may  be  easily  calculated  for  various  host  environments,  a  modi¬ 
fied  EAM  may  be  used  to  analyze  the  impurity’s  energetics  in  a  variety  of  situations. 
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